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UNLIMilED 


LIST  OF  SYMBOLS 


Cartesian  coordinates 
Cylindrical  polar  coordinates 
Spherical  polar  coordinates 
Three-dimensional  Green's  function 
Three-dimensional  free-space  Green's  function 

=G_(x,y,z,x  ,y  ,z  )  in  Cartesian  coordinates 
r  O  O  O 

=G„(r,4>,z,r  ,<P  ,z  )  in  cylindrical  polar  coordinates 

r  O  O  O 

=Gp(R,9,<t),R^,0^,‘j>^)  in  spherical  polar  coordinates 

Two-dimensional  free-space  Green's  function 

=G„(x,y,x  ,y  )  in  Cartesian  coordinates 
t  o  o 

=Gp(r,^,r^,(>^)  in  polar  coordinatinates 

One-dimensional  free-space  Green's  function 
Three-dimensional  Dirac  delta  function 
One-dimensional  Dirac  delata  function 
Bessel  function  of  the  first  kind,  order  n 

Hankel  function  of  the  first  kind,  order  n,  H  (x)=J  (x)+iY  (x) 
Bessel  function  of  imaginary  argument,  order  n 
Spherical  Bessel  function  of  the  first  kind,  order  n 
Note  j^(x)='/'(Tr/2x)J^^l^^  Cx) 

Legendre  function  of  the  first  kind,  degree  n 

Associated  Legendre  function  of  the  first  kind  of  degree  n 

and  order  m.  Note  P™  ■P°  /  P™=0,  m>n 

n  -(n+1)  n 

e^=l;  e^=2,  n=l,2,3,  ... 

(-u:/c)  wavenumber 
»(z-z^)sgn(z-z^) 

[(x-X^)2+(y-y^)2+(z-z^)2] 

=/  Ir^+r^  -  2rr  cosCiji-^!  )  +  (z-z  )^] 

O  O  O  o 

[  r2+r2  -2R^Rcos(ij;)l  ,  cos(i|;)»  cosQ  cos9^  +  sin0  sin9^  cos (<|)-<))^) 

[  (x-x^)2+(y-y^)2] 


W  {  r2+r^-2rr  cosC'J-'t’  )) 


INTRODUCTION 


The  Green's  functions  of  the  reduced  wave  equation  are  defined  as  the 
solutions  of  the  inhomogeneous  equation 

(v2+k2)G(R,R^)  =  -4it6(R-R^)  (1.1) 

which  satisfy  the  appropriate  linear  boundary  conditions.  In  the  absence  of 

boundaries  the  solution  is  the  free-space  Green's  function,  Gp,(R,R  ),  which 
takes  the  particularly  simple  form  ° 

Gj.(R,R^)  =  exp(ik|R-R^j)/lR-R^|  (1.2) 

for  outgoing  waves  in  three-dimensions.  The  time  dependence  factor,  exp(-i(i)t), 
which  is  assumed  throughout,  is  omitted  from  all  equations.  It  is  evident 
that  this  Green's  function  is  the  outgoing  wave  solution  of  the  reduced  wave 
equation  when  the  excitation  is  a  point  source.  It  is  also  evident  that 
reciprocity  applies,  viz  G^  is  unaltered  when  the  positions  of  the  source  and 
the  observer  are  interchanged. 

The  reduced  wave-equation  with  general  source  distribution  F(R)  inside  a 
volume  V 


(72+k2)p(R)  -  -4ttF(P  (1.3) 

is  easily  solved  by  the  principle  of  superposition  when  the  Green's  function  is 
known.  The  solution  is 

p(p  -  I  F(R')G(R,R')d3R'  (1.4) 

J  V» 

The  determination  of  a  Green's  function  usually  proceeds  by  assuming  that 
it  is  the  sum  of  the  free-space  Green's  function,  G^.,  and  a  scattering  term 
which  satisfies  the  homogeneous  equation  and  whose  amplitude  is  determined  by 
the  boundary  conditions.  In  this  case  the  homogeneous  form  of  equation  (1.1) 
is  usually  solved  by  the  method  of  separation  of  variables,  which  results  in 
the  scattering  term  being  represented  by  integral/series  transforms.  There  is 
therefore  a  requirement  to  expand  the  free-space  Green's  function,  G  ,  in  a 
i  similar  form  in  order  to  find  the  amplitud  of  the  scattered  wave. 


The  expansions  of  the  free-space  Green's  function  in  the  various  coordinate 
systems  are  developed  in  texts  too  numerous  to  reference  here.  It  is  the 
purpose  of  this  memorandum  to  obtain  the  expansions  in  a  simple  way,  without 
being  especially  rigorous.  It  is  hoped  that  the  expansions  will  be  of 
particular  interest  to  researchers  who  are  using  Green's  function  techniques 
for  perhaps  the  first  time,  or  to  those  who  would  find  a  collection  of  free- 
space  Green's  function  representations  to  be  useful.  For  reference,  the 
Appendix  summarizes  the  representations  derived  herein. 
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2.  EROBLEM  FORMULATION 

The  free-space  Green's  function  considered  here  is  the  outgoing  wave 
solution  of  equation  (1.1),  which  is  expressed  in  Cartesian,  cylindrical 
and  spherical  coordinate  systems  respectively.  Figure  1,  as 

[  92/ 3x2+32/ 3y2+32/92'2+ij2]G(R  R  )=-4it6(x-x  )6(y-y  )6(z-z  ) 

“  — o  o  o  o 

(  92/  3r2+  (1/r)  9/  9r+  (l/r2)  92/  9*2+32/ 222+.k2]  G.(R,^) 

*  -An6(r-r^)S(<t>''4>^)6(z^z^)lr 

[  92/3R2+(2/R)9/9R+(l/R2)92/9e2+(cot  e/R2)3/9e 

+  (l/R2sin26)92/3<|>2+k2]G(R,R  ) 

—  — o 

=  -4ir6  (R-R^)  6  (9-0^)  6  (<1)-(|)^)/R2sine 

The  evaluation  of  the  three-dimensional  free-space  Green's  function  in 
each  coordinate  system  proceeds  as  follows.  First,  the  differential 
equation  (1.1)  is  solved  by  the  method  of  separation  of  variables 
which  leads  to  a  general  integral/series  transform  representation  of 
the  Green's  function.  Secondly,  this  transform  representation  is 
evaluated  to  obtain  the  Green  i  function  in  a  form  appropriate  to 
outgoing  waves .  Thirdly,  the  transform  representation  is  evaluated  to 
recover  the  Green's  function  in  its  simple  form,  equation  (1.2). 

Finally  the  particular  cases  of  the  two-dimensional  and  one-dimensional 
free-space  Green's  function  are  analysed. 


3.  CARTESIAN  COORDINATES 


3.1  Fourier  Transform  Solution 


The  three-dimensional  Fourier  integral  transform  representation 


Gp(R,R^)  -  (1/8it^)  j  Gp(a,6,e,R^)exp  (  iox+iSy+iezJ  dotdSde 


(3.1) 


and  its  inverse 


Gp(a,B,e,  R^  -  j  Gp(R,R^)exp  [ -iax-i^-riexj  dxdydz 


(3.2) 


are  used  to  obtain  a  general  solution  of  the  reduced  wave  equation  (2.1), 
with  point  source  excitation,  in  Cartesian  coordinates  as 


I  exp  (  ia(x-x  )  +  ig(y-y  )  +  ie(z-z  )]dod6de  (3.3) 

- 2 - 2 - 2 - 

a2+62+e2.k2 


The  Green's  function  and  its  Fourier  transform  G^  are  considered  as 
distributions  in  the  sense  of  Schwartz  [4]  because  consideration  must  be 
given  to  the  choice  of  integration  contours  in  their  evaluation  to  ensure 
that  the  outgoing  wave  solution  alone  is  obtained. 


3.2  Plane  Boundary,  z=constant 

In  many  cases  of  interest  there  are  plane  boundaries,  perpendicular  to 
the  z-axis,  say,  on  which  boundary  conditions  must  be  satisfied.  In  such 
cases  it  is  advantageous  to  carry  out  the  c-integration.  This  integration 
may  be  accomplished  by  choosing  a  contour  which  depends  on  the  sign  of 
z-z^,  viz  for  z-z  ^0  the  contour  is  closed  in  the  upper  half  plane,  while 
for  z-z^  <  0  the  contour  is  closed  in  the  lower  half  plane.  These 
contours  are  shown  in  Figure  2. 


When  k^>  +  3^  the  integrand  has  poles  on  the  real  axis  and  the  portion 

of  the  contour  along  the  real  axis  must  be  indented  to  avoid  passing 
through  them.  The  integral  representation  of  the  Green's  function, 
equation  (3.3),  is  a  general  solution  which  is  arbitrary  with  respect  to 
a  solution  of  the  homogeneous  wave  equation  corresponding  to  incoming 
free-waves  generated  at  infinity.  It  is  necessary  to  indent  the  contour 
so  that  outgoing  waves  only  are  present.  This  'radiation  condition'  may  be 
simulated  mathematically  by  allowing  for  energy  dissipation  in  the  medium 
by  setting  k=k+ik',  corresponding  to  a  complex  sound  velocity  c=c-ic',  so 
that  any  disturbance  originating  at  infinity  is  negligible  in  any  finite 
region  of  interest.  The  'real'  poles  are  thus  forced  off  the  real  axis  to 
allow  an  unambiguous  evaluation  of  the  integral  by  the  residue  theorem. 

The  dissipation  term  is  then  allowed  to  tend  to  zero  to  show  that  this 
k '-prescription  is  equivalent  to  indenting  the  contour  above  and  below  the 
negative  and  positive  poles  respectively.  The  following  integrals  are 
obtained 

00  90 


Gj.(R,R^)  =  (i/2Tr) 


J 


(l/Y)exp  [  ia(x-x^)  +  igCy-y^^)  +  iY(z-z^)]  dadg 


—  00—00 
00  00 


2  >  *  (3.4) 

o 


(1/Y)exp  [  ia(x-x^)  +  iS(y-yjj)  -  iY(z-z^)]  dadB 


—00  —00 


Z  C.  X 


(3.5) 


or  more  generally 

Gp(R,R  )  -  (i/2ir) 
r  "O 


*00  —00 


(1/Y)exp  (  ia(x-x^)  +  i8(y-y^)  +  iY!z-eQ|] dodg 

(3.6) 


for  all  z,  where 

Y  ■  ,  Im(Y)^»  Re(Y)>o 
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3 . 3  Reduction  to  Standard  Form 


The  standard  form  of  the  Green's  function  equation  (1.2)  may  be 
obtained  by  evaluating  the  double  integral  in  equation  (3.6). 

This  is  accomplished  by  first  transforming  the  variables  to  polar 
coordinate  systems,  via  the  substitutions 

a  “  s.cos  6  3  ••  s.sin  9  dad6  *  sds  d0 

x-x  ■  t.cos  <|)  y-y  ■  t.sin  it 

The  identity  [2] 


[  its.cos(6-(ji)l  de  «2irJQ(ts) 


is  then  used  to  reduce  the  resulting  double  integral  to  the  single 
integral 


CO 

/ 


J  (ts)exp  t -Iz-z  |N/’(s^-k^)]  sds 
o _  o 

^ (s^-k^) 


(3.7) 


The  integrand  has  a  branch  point  at  s*k,  and  the  contour  must  be 
indented  to  avoid  it.  The  inclusion  of  a  dissipation  term,  as  in 
the  case  with  a  plane  boundary,  allows  the  branch  point  to  be  lifted  above 
the  real  axis.  In  the  limiting  case  of  no  dissipation  this  is  equivalent 
to  indenting  the  contour  below  the  branch  point.  Figure  3.  With  this 
choice  of  contour  for  the  integration,  ensuring  that  the  outgoing  wave 
solution  alone  is  obtained,  equation  (3.7)  is  a  standard  integral  [1] 
which  reduces  to 


exp  [ik{(x-x  )2  +  (y-y  )2  +  (z-z 

G  (R*R  )  - - 2 - 2 - SL - 

(  (x-x^)2  +  (y-y^)2  +  (z-z^)2] * 

which  is  the  simple  form  of  the  free-space  Green's  function  in  Cartesian 
coordinates 

4.  CYLINDRICAL  COORDINATES 
4.1  Transform  Solution 


The  Green's  function  is  represented  by  a  Fourier  series  transform  in 
the  angle  (p,  a  Fourier  integral  transform  in  the  axial  coordinate  z, 
and  a  Hankel  transform  in  the  radial  coordinate  r,  viz 


Gj,(R,R^)  -  (l/2Tr)Z  exp(in((»)  | 


n*“oo 


G_(C,n,a,R^)J  (5r)exp(ioz)5d5da 
^  “O  n 


j  I 
—00  0 


(4.1) 
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where  the  inverse  transform  is 

Pf»  2^  00 


Gp(5,n.ci,R^)  *  (l/2ir) 


Gp(R,R^)J^(5r)exp(-in^-iaz)rdrd<tidz  (4.2) 


-00  0  0 


The  above  transform  representations  are  used  to  obtain  a  general  solution 
of  the  reduced  wave  equation  (2.2),  with  point  source  excitation.  After 
simplification,  using  the  differential  equation  satisfied  by  the  Bessel 
function,  this  general  solution  is 


n=— 00 


; 

o 


J^(5r^)J^(Cr)exp  [  ia(z-z^)] Cd^da 


(4.3) 


4.2  Cylindrical  Boundary  r~constant 

In  many  cases  of  interest  in  this  geometry  there  are  cylindrical 
boundaries  present,  defined  by  r=constant,  on  which  boundary  conditions 
must  be  satisfied.  In  such  cases  it  is  advantageous  to  carry  out  the 
^'integration.  The  integrand  has  a  pole  on  the  real  axis  for  k^>a^, 
and  the  contour  must  be  indented  to  avoid  it.  Introducing  a  dissipation 
factor,  as  before,  to  obtain  the  outgoing  wave  solution  only,  allows  the 
pole  to  be  lifted  above  the  real  axis  or  equivalently  the  contour  to  be 
indented  below  the  pole.  The  standard  integral  [1] 


► 

■’o 

with  Im(Y*v'^(k^'Oi^))>0, 
represented  as 


(Tri/2)J^(Yr^)H^(Yr)  r  >  r^j>0 

*  < 

(Tri/2)Jn(Yr)H„(Yrjj)  r^>  r  >0 

then  enables  the  free-space  Green's  function  to  be 


Gp(R.Ro) 


(i/2)  I  exp  I  in(<ti-(|i  )] 
n»-®  ° 


(i/2)  Z  exp  (  in(<Ji-<( 


Qm— 00 


00 

'"I 


J  (Yr„)H  (Yr)exp  [  ia(z-z  )]  da 
non  o 

r  >  r  >0  (4.4) 

o 

Jjj(Yr)Hj^(Yr^)exp  [  ia(z-z^)]da 
r  >'T  >0 


with  Y■^/'(k2-  a^)  and  with  Im(Y)>0,  Re(Y)^. 


(4.5) 


4.3  Source  or  Observer  on  z-axis 


In  the  special  case  of  the  source  being  positioned  on  the  z-axis,  it  is 
evident  that  the  Green's  function  must  be  axisymmetric.  Hence,  the  n^O 
term  alone  contributes  to  the  field  which  is  obtained  from  equation  (4.3) 
as 


Sr(5.S„) 


J^(Cr)exp  [  ia(z-z^)]  Sd^do 


(4.6) 


The  ^-integration  is  facilitated  by  extending  the  range  of  integration 
to  the  entire  real  axis  by  the  use  of  the  formula 

J  (Cr)  =  (1/2)  [H^(5r)-H^(-Cr)]  (4.7) 

o  o  o 


to  give 


G  (R,R  )  »  (l/2ir) 


— flO  —00 


H^(Cr)exp  [  io(z-^)} ^d^da 


(4.8) 


Figure  3  shows  the  contour  used  to  evaluate  the  ^-integral.  The 
indentation  at  the  origin  is  necessary  because  of  the  branch  point  of 
the  function  H  (Cr) .  The  contribution  from  this  indentation  can  easily 
be  shown  to  be”zero  when  the  small  argument  expansion  of  the  Hankel 
function,  viz  Uo(x)=>(2i/ir)  jln(x) ,  is  used  in  its  evaluation.  The 
integrand  also  has  poles  on  the  real  axis  when  The  required 

indentations  are  again  chosen  by  the  inclusion  of  a  small  dissipation 
term,  k=k+ik^,  to  give,  on  applying  the  residue  theorem. 


00 


Gp(R,R^)=  (i/2) 


H^(Yr)exp  [  ia(Z“Z^)] da 


(4.9) 


a  result  which  is  identical  to  equation  (4.4),  in  the  special  case  of  >^^=0. 

In  the  case  of  the  observer  being  positioned  on  the  z-axis  the  analysis 
is  identical,  as  would  be  expected  by  reciprocity,  and  yields 


(i/2) 


H  (yr  )exp  [  ia(z“z  )]do 
o  o  o 


(4.10) 


Equations  (4.4)  and  (4.5)  are  therefore  valid  for  all  values  of  r  and  r 
such  that  r  and  r  are  not  both  zero. 


4 . 4  Plane  Boundary  z=cons tant 


In  cases  where  a  plane  boundary  perpendicular  Co  Che  z-axis  is  presenc, 
on  which  boundary  condicions  muse  be  sacisfied,  ic  is  advancageous  Co 
perform  Che  a-inCegracion  in  equacion  (4.3).  The  concours  used  for 
chis  are  idencical  Co  chose  used  for  Che  e-incegracion  in  Seccion  3, 
and  give  Che  resulc 


00 


®  i  2  exp  [  in((j)-({i  >] 

r  o  Q 

na-oo 


) 

o 


Jri(5r)J^(Cro)exp  [  i [ z-z^ | Cd5 


(4.11) 


4.5  Reduccion  Co  Scandard  Form 

Equacions  (4.4)  and  (4.5)  can  be  reduced  Co  sCandard  form,  equacion  (1.2), 
by  making  use  of  one  of  Che  Bessel  funecion  addicion  cheorems  [3]  ,  viz. 


H  (Yri) 


Z  exp  [  in(iji-!j>  )]  H  (yr  )J  (yr) 
n=-«o  o  n  o  n 

y  exp  [  in(<|i-<j)  )]  H  (yr)J  (yr  ) 
L  o  n  n  o 

n=-«o 


Co  give 


r  >  r>0 
o 

r  >  r  >0 
o 


Ho(Yri )exp  [  ia(z-z^)] do 


(4.12) 


where 

r^  »  >/■  [  r^  +  r|  -  2r  r^  cos  (i)>-())^)]  *  ■J  [  (x-x^)^  +  (y-y^)^] 

This  scandard  inCegral  [  2,p.710j  reduces  Co  Che  simple  form  of  Che 
free-space  Green’s  funeCion  for  ouCgoing  waves,  equacion  (1.2). 


SPHERICAL  COORDINATES 


5.1  Transform  SoluCion 


The  Green's  funecion  is  represenced  by  a  Fourier  Series  cransform  in  Che 
angle  <t)  ,  an  associaCed  Legendre  Series  Cransform  in  Che  angle  6,  and  a 
spherical  Hankel  cransform  in  Che  radial  coordinace  r,  viz 


00  00 
■  Z  exp(imo)  Z 
ni=  -00  n=— « 


?®(cos  3) 
n 


r 


G  (5,n,m,R  )j  (?R)C^/2^C 
r  On 


where  the  Inverse  transform  is 

G  (C,n,m,R  )  =  (2n+l)  (n-m)  I)  . 

r  — o 

2ir  ir  « 

*  f  f 

G„(R,R  )exp(-im<|i)j  (5R)P“(cos  e)R2sin0  dRd0d(J»  (5.2) 

P - o  n  n 

) 

0  0  0 

The  above  transform  representations  are  used  to  obtain  a  general  solution 
of  the  reduced  wave  equation  (2.3)  in  spherical  coordinates  which,  after 
simplification  using  the  differential  equation  satisfied  by  the  spherical 
Bessel  function,  is 


00 

(i/tt)z: 

I^s^oo 


ln|  +  1 
Z 

m=-  I  n  j  -.1 


(2n+l)  (n-m)!  p“(cos  0)p“(cos  0Q)exp  [im((|i-(()Q)] 
(n+m)! 


00 


0 


j  (CR  )j  (5R)c2d5 
n  On 


(5.3) 


5.2  Spherical  Boundary 

In  cases  of  interest  in  this  geometry  there  are  spherical  boundaries 
present,  defined  by  R=constant,  on  which  boundary  conditions  must  be 
satisfied.  In  such  cases  it  is  advantageous  to  carry  out  the  ^'integration 
in  equation  (5.3).  The  identity 


'o 


2v/’(RR  ) 
o 


C2-k2 


(5.4) 


enables  this  integration  to  be  accomplished  by  the  methods  of  Section  4  to 
give  the  outgoing  wave  representation  for  0<R^<R 


Gp(R.Ro)  >  [iri/4>^(RR^)]Z  (2n+l)  ^kR)  . 

n“— oo 


m* 


inl-Hl 

Z 

-|n|-l 


(n-m)lP™(cos  9)P™(cos  0^)exp  {  im((J)-(J)^)] 
(n+m) I 


(5.5) 


while  for  0<R<R  it  is  only  necessary  to  interchange  R  and  R^.  The  above 
equation  may  be  simplified  to 


Gp(R.R^)  =  [Tri/2/(RR^)]  Z  (2n+l) ,  (kR)H^^l .  (kR^) . 

n=o  ^ 

n  (n-m)!  P  ™  (cos  0)P™(cos  0  )  cos  [in(4>-<i>  )]  (5.6) 

j,  _ n _ n _ o _ o 

m=o  (n+m)! 


5.3  Source  or  Observer  at  the  Origin 

In  the  special  case  of  the  source  being  positioned  at  the  origin  it  is 
evident  that  the  Green's  function  is  spherically  symmetric.  Hence  the 
(n=0,  in=0)  and  (n=-l,  m=0)  terms  alone  contribute  to  the  field  which  is 
obtained  from  equation  (5.3)  as 


Gp(R,Ro) 


(2/Tr) 


j^(4R)c2d5 

5^  -  k2 


(5.7) 


This  integral  is  of  standard  form  when  the  spherical  Bessel  function  is 
replaced  by  a  Bessel  function  of  order  one-half  and  the  contour  is 
indented  below  the  pole  at  C=h.  It  is  [2,  p.687] 


Gp(R,R^)  =  /(-2ik/TrR)Ki^j,(-ikR) 
which  can  be  expressed  [2,  p.967]  as 

G„(R,R  )  =  exp(ikR)/R  (5.8) 

r  • 

In  the  separate  case  of  the  observer  being  positioned  at  the  origin,  the 
analysis  is  identical  and  gives  the  result  expected  by  reciprocity,  viz 
equation  (5.8)  with  R  and  R^  interchanged. 

5.4  Reduction  to  Standard  Form 

Equation  (5.6)  can  be  reduced  to  the  standard  form  of  the  Green's  function 
equation  (1.2)  by  making  use  of  the  addition  theorems  [2,  p.l013  and  p.980] 
for  the  Legendre  and  spherical  Bessel  functions,  viz 


n  e  (n-m)!  P  (cos  0)P™(cos  0  )cos  [m(d)-(j)  )]  (5.9) 

P^(cos  ip)  =  Z  -H - II - ? - 2 - 2_ 

m=o  (n+m)! 


exp(ik|R-R^|)/lR-R^l  =  [  ui/2/(RR^)]  Z 

R  >  R  (5.10) 


n=o 


where  cos  ip  =  cos  0  cos  0  +  sin  0  sin  0  cos  ) 

0  o  o 


is  the  cosine  of  the  angle  between  the  source  and  observer. 


ONE  AND  TWO  DIMENSIONAL  REFRESRNTAtTONS 


6.1  General 

The  analysis  necessary  to  obtain  the  expansions  of  the  one  and  two 
dimensional  free-space  Green's  functions,  and  their  subsequent  reduction 
to  the  simple  standard  forms  follows  the  same  procedure  as  that  used  to 
solve  the  three-dimensional  problem.  This  section  therefore  contains  a 
minimum  of  analytical  detail. 

6.2  Cartesian  Coordinates;  (x.y) 

The  two-dimensional  Fourier  transform  and  its  inverse  are  used  to  obtain 
the  general  representation  of  the  free-space  Green's  function  in  two- 
dimensions  as 


00 


00 


Gp(r,r  )  -(l/ir) 
r  —  — o 


exp  (  ia(x-x^)+i6(y-y^)]  dode 

o2+g2_jt2 


(6.: 


In  the  presence  of  the  plane  boundary  y-constant,  the  3-integration  is 
performed  to  give  the  free-space  Green's  function  as 


<9 


exp  [  ia(x-x^)  +  i^Ck^-a^) Jy-y^ )]  da  (6. 


This  integral  may  be  written  as 


cos  { a(x-x^)lexp(-|y-y^|/(a2  +  (-ik)2)ldo 
>A(a^  +  (-ik)^) 

0 


which  is  of  standard  form  [2,  p.498]  when  damping  in  the  system  is 
considered,  and  gives 


Gp(r,r^)  -  2K^(-ik  I  (x-x^)2  ♦(y-y^,)^)  ^2) 

which  reduces  to  the  free-space  Green's  function  in  its  simple  form,  viz 

GpClff  )  ■  niH  (k  [  (X-X  )2+(y-y  )2]  1/2)  (g  ' 

r  o  o  o  o 

-  iriH^(k|r-r^|)  (g., 


-  14  - 


6.3  Polar  Coordinates;  (r,(|)) 

The  Fourier  series  and  Hankel  transforms  in  the  angle  ({>  and  the  radius 
r,  respectively,  and  their  inverses  are  used  to  obtain  the  general 
representation  of  the  free-space  Green's  function  in  two  dimensions  as 


)  =21  exp  [  in((t)-(j)  )] 
r  o  o 


na-oo 


J  (?r)J  (Cr  )5dC 
n  n  o 


(6.5) 


52-k2 


In  the  presence  of  a  circular  boundary,  tne  ^'integration  is  performed  as 


Gp(r,r^)  *  Tri  I  exp  (  in ((ti-<Ji^)]  H  (kr)J  (kr^) 
n»-<»  ^ 

00 

Gp(r,r  )  =  iri  Z  exp  [  in(ij)-<t)  )J  H  (kr  )J  (kr) 
t  o  o  n  o  n 

n=-« 


r  >  r  >  0  r+0 

o 

(6.6) 

r  >  r  >0  r  +0 
o  o 

(6.7) 


a  result  which  reduces  immediately,  via  a  Bessel  function  summation  theorem 
[  3]  ,  to  the  simple  form  of  the  two-dimensional  free-space  Green's  function 
equation  (6.4). 

6.4  One-Dimensional  Representation 

The  one-dimensional  Fourier  transform  and  its  inverse  are  used  to  obtain 
the  general  representation  of  the  free-space  Green's  function  in 
one-dimension  as 


Gp(x,Xo) 


exp  I  iQi(x-x^)]  /  (a2  -  k2)dri 


(6.8) 


which  is  immediately  integrable  to  the  standard  and  simple  form  of  the  one¬ 
dimensional  free-space  Green's  function,  viz 


Gp(x,x^)  =  (2Tri/k)exp  (  ik(x-x^)l  x  >  Xq  (6.9) 

Gp(x,x^)  “  (2Tri/k)exp  [  ik(xo-x)J  x  <  Xq  (6.10) 


E  A  Skelton  (SO) 
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APPENDIX 


GREEN'S  FUNCTIONS  IN  THREE  DIMENSIONS 
General  Form 


“  exp(ik|R-R^|)/|R-R^ 


GpCR.Ro)  »  (l/27r2) 


exp[  ia(x-x  )  +  iS(y-y  )  +  ie(z-z  )ldad6de 
_ o _  o _ o _ _ 

a2+g2+e2_ij^2 


G„(R,R  )  =  (I/tt)  Z  exp  [  in((|)-(|i  )] 

r  “  “O  O 

n=-<» 


J^(5r)J^(5r^)exp  [  ia(ZT.z^)l  Cd^da 


n  1+1 


G„(R,R  )  -  (I/tt)  Z  I  ,  ^ 

'  ®  n=-«o  m=-ln|*-l 


(2h+l)  (n-m)!  P™  (cos  0)P°  (cos  6  )exp  [  im( 4>-i^ 
n  n  o  o 


Cartesian  Coordinates  -  Plane  Boundary  (z=constant) 
Gp(R,R^)  -  (i/Zir) 


(1/Y)exp  [ia(x-x^)  +  iS(y-yjj)  +  iyC  z-z^)]  dadS 


with  y  “  '^(k2-a2-g2j  ^  £qj.  (^-z  )>0.  Interchange  z  and  z  when  (z-z  )<0. 


Cylindrical  Coordinates  -  Cylindrical  Boundary  (r«constant) 


Gp(R,R^)  -  ^ 


(i/2)  Z  exp  [  in((Ji-(j)  )J 
n-  -«o  0 


J  (yr  )H  (Yr)exp  [  ici(z-z„)]do 
non  o 


*00 

00 

r 


(i/2)  Z  e  cos[  n(^-<^  )] 


J  (yr  )H  (Yr)exp[  ia(z-z  )]  do 


)]  /Cn+m) 


i  L  exp[ 
n=-“ 


J  (5t‘)J  (5r-)exp  [  i(z-z  )'^(k^-5^)l 
n  _ n  o  o _ 


/(k2-c2) 


i  Z  e„  cosl  n(((i-(j  )] 
n  o 

n®o 


J^(5r)J^(Cr^)exp  [iCz-z^)'/’ (k^-o^)]  5d5 


for  (z-z  )^.  Interchange  z  and  z  when  (z-z  )<0. 
o  o  o 


Spherical  coordinates  -  Spherical  Boundary  (R=constant) 


n|  +  l 


Z  (n-m)!  p”  (cos  8)p”  (cos  0^)exp  {  im(i^-4  )J /(n+m)! 

M.  n  no  0 

-1 


<=r<S-So)  -i 


h”0 


E  (n-m)!  P®  (cos  9)P®  (cos  9  )cos  [in((^-^  )]/(n+m)! 
n  n  o  o 

ni“o 


for  0<R<R^.  Interchange  R  and  R^  when  0<R^<R. 


GREEN '.S  FUNCTIONS  IN  TWO  DIMENSIONS 


General  Form 


S<E'Eo>  ■ 


exp  [  ia(x-x^)  +  i0(y-y^)] dadg 


a2+62_k2 


®F^E»Eo)  “21  exp[  in((j)-(|i  )] 


n" 


J„(er)J„(«r  ){dC 


r2_i,2 


"it  "J  %  -.  ■:^’ L"»  V  * 


Cartesian  Coordinates  -  Plane  Boundary  (y=constant) 


Gp(r.r^)  »  i 


exp  (  ia(x-x^)  +  i(y-y^)^/(k^-ol^)l  da 
/ (k^-a^) 


for  y^Q*  Interchange  y  and  y^  when  y  ^y^- 


Polar  Coordinates  ~  Circular  Boundary  (r=constant) 
f 


iri  Z  exp  [  in(<|)“4i^)]  H^(kr)  J^(kr^) 
n«  —00 


Tri  Z  e  cos  (nCdi-A  )1  H  (kr)J  (kr  ) 
n=o  ^  n  no 


for  oTq*  Interchange  r  and  r^  when 


GREEN'S  FUNCTION  IN  ONE  DIMENSION 


General  Form 


Gj,(x,x^)  -  (2Tri/k)exp(ik|x-x^| ) 


exp  I  ia(x-x^)]  da 


-  21 


!<*>•<*+  8* 


Z  -  Zo  ^  0 


END 
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